Temperature dependence of the conductivity of ballistic graphene 
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We investigate the temperature dependence of the conductivity in ballistic graphene using Lan- 
dauer transport theory. We obtain results which are qualitatively in agreement with many features 
recently observed in transport measurements on high mobility suspended graphene. The conduc- 
tivity a at high temperature T and low density n grows linearly with T, while at high n we find 
a ~ \/\n\ with negative corrections at small T due to the T-dependence of the chemical poten- 
tial. At moderate densities the conductivity is a non-monotonic function of T and n, exhibiting a 
minimum at T = 0.693hvy^]n\ where v is the Fermi velocity. We discuss two kinds of Fabry-Perot 
oscillations in short nanoribbons and their stability at finite temperatures. 

PACS numbers: 72.10.-d,73.50.-h 



Ballistic transport in graphene has been theoreticaUy an- 
alyzed P, 13, ^ shortly after the pioneering measure- 
ments of the quantum Hall effect in the first single- 
layer graphene samples 0, @. So far theoretical ap- 
proaches concentrated on low temperatures and density, 
predicting a finite and universal minimum conductiv- 
ity a = (4/7r)e^//i as well as a universal Fano factor 
F = S/2eI —1/3 (the ratio of the current noise S and 
the average current /, e > being the electron charge). 
This is a surprising result in the ballistic transport regime 
where one might naively expect no current noise at all. 
The reason for finite noise lies in the dynamics of charge 
carriers in graphene which at low energies is governed 
by the massless Dirac equation and not by the massive 
Schrodinger equation. As a result, the ballistic transport 
resembles diffusive transport in a normal metal, which 
has been coined pseudo-diffusive transport in the litera- 
ture. Both the predictions for the conductivity as well 
as the Fano factor have been experimentally observed 

Recently, mobilities approaching 200000 cm^/Vs have 
been reported for ultraclean suspended graphene 0, [l3| 
which are an order of magnitude larger than typical mo- 
bilities of graphene deposited on a substrate. The trans- 
port characteristics of these experiments suggest that the 
samples reach the ballistic regime with respect to disor- 
der scattering. Indeed, the conductance scales with the 
number of channels, which is inconsistent with dominant 
scattering from charged impurities, or ripples llj, the 



latter being the most likely elastic scatterers in suspended 
graphene. Further, a is proportional to the sample length 
L for small sizes and low temperature T. 

While these two features suggest balHstic transport, they 
do not actually rule out the presence of inelastic scatter- 
ing due to electron-electron interactions, which preserve 
both the density dependence a ~ -^/n and the propor- 
tionality to L [12] - except for n — and finite T where 
interactions lead to a L-independcnt finite conductivity. 



in contrast to the ballistic case (see Eq. ([7]) below). The 
data of Ref. Q is indeed likely to bear fingerprints of 
Coulomb interactions in a certain parameter regime close 
to the neutrality point and at higher temperature where 
interactions are strongest |13j . On the other hand, how- 
ever, the significant linear increase of <t{T) at charge neu- 
trality reported in Ref. [lol| is inconsistent with dominant 
electron-electron interactions. (In the experiment [loj 
the latter are presumably screened by nearby metallic 
electrodes.) This insulating trend in the temperature de- 
pendence, and the opposite metallic trend at moderate 
density, which is reported in both Refs. 0, , has re- 
mained a puzzle. The authors of Ref. P attempted to 
explain the latter by scattering from phonons. However, 
they point out that the decrease of the effect with den- 
sity is not consistent with such a scenario. Ascribing 
the observed effects to electron-electron interactions is 
not consistent either, since those would exhibit a rather 



weak temperature dependence 13|, [lj|, with a metallic 
instead of an insulating trend at low density. Hwang and 
Das Sarma [3] have proposed scattering from charged 
impurities and their temperature dependent screening as 
a possible explanation for the different behavior at low 
and high density. However, this model cannot account 
for the density dependence of the data and the length 
independence of the conductance at low temperature. 

In this Letter, we analyze impurity-free transport at fi- 
nite T. As motivated above, we assume interactions to 
be weak, either due to dielectric attenuation from a sub- 
strate, or screening by nearby metals or the finite density 
of electrons themselves. We show that the extension of 
the ballistic transport model of Ref. to T > qualita- 
tively explains most of the features observed in Ref. [l3| 
and, apart from the low n and high T regime, also those 
in Ref. In particular, in the absence of interactions, 
the minimum conductivity is shown to grow linearly with 
T. At finite carrier density n, the initial T dependence 
is negative, but changes sign at a temperature of order 
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T ~ hv^/\n\ where v « 10^ m/s is the Fermi velocity. 
In the model for balhstic graphene proposed in Ref. [1], 
a nanoribbon of width W is suspended between left and 
right reservoirs (wide graphene regions), which are a dis- 
tance L apart. The leads and the sample region are sub- 
ject to a step-like electrostatic potential 



0(x) 



a; < (left lead), 
< X < L (sample), 
X > L (right lead). 



(1) 



the zero of energy being chosen as the Dirac point in 
the sample region. Modeling the leads by highly doped 
graphene {\4>oo \ much larger than the gate induced Fermi 
level in the sample) one finds that the parameter 
drops out in the end. 

Ballistic transport requires the sample length to be 
shorter than the mean free path £. In nearly impu- 
rity free samples the latter is limited by inelastic scat- 
tering, dominated at low T by Coulomb scattering, 
£ « (fi?;/a2)T2/max(|/i|,r) [ij, Q. The latter is 
strongest in the non-degenerate regime T > but is 
relatively weak outside because of screening effects, and 
the down-renormalization of the Coulomb coupling con- 
stant a. The latter is of order 0(1) in the unscreened 
situation pertaining to free-hanging graphene in the non- 
degenerate regime. However, on a clean substrate with 
large dielectric constant, or in the presence of nearby 
metallic contacts, as in Ref. the effective value of a is 
substantially reduced. This opens a reasonably large win- 
dow of applicability for the non-interacting theory which 
we develop below. 

The zero temperature conductivity at a fixed Fermi level 
E in the sample is given by the Landauer formula 



ME) = ^GiE) 



W h 



(2) 



n=0 



where G{E) is the conductance, g = A is the degener- 
acy due to spin and valley degrees of freedom, and n la- 
bels the transverse modes of the graphene ribbon. Their 
transmission probability follows from solving the propa- 
gation through the potential ([T|), 



Tn{E) 



E"^ — [hvqnY cos2(fc„L) ' 



(3) 



where fc„ = [hv)~^ E"^ — (hvqn)^. The transverse mo- 
mentum q„ is defined for various boundary conditions as 
qn = [n + 7)7r/T/F. Below we use 7 = 1/2 corresponding 
to infinite mass confinement, see Ref. [2]. 
The energies Ew,l = hv/{W, L} set typical scales below 
which finite size effects are important. At higher energies 
confinement effects are unimportant. This is easily seen 
in the asymptotics of the T = conductivity in the limit 



W ^ L (which is independent of 7 Q) 
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I [l-l-0.10094?/2 + O(y4 



y < 1 



(4) 



y + 



sin(2j^ — 7r/4) 



+ 0(l/y), 



This reproduces the well-known minimal conductivity 
CTmin(M = r = 0) = ^i. For E > El, we find a 
to be essentially proportional to L and E, which refiects 
the number of conducting channels at energy E. The 
oscillatory interference term will be discussed in detail 
further below. 

We now take into account the effects of finite temperature 
in this model. In linear response, the Landauer formalism 
yields the exact formula for the conductivity 

a{^l,T)=T-' I ao{E)f{E)[l^f{E)]dE, (5) 

where f{E) = (1+cxp [{E - ^)/r])"i is the Fermi distri- 
bution function (fcs = 1 throughout) and /i is the chem- 
ical potential in the sample. 

For T ^ El it is justified to neglect the oscillatory term 
in ([4]), and we find the T-dependent conductivity 



h hv 



|Ai|+2riog(l 



) • (6) 



At charge neutrality, and in the non-degenerate case, 
T ^ the conductivi ty g rows linearly with T, simi- 
larly as observed in Ref. 1C||, 



LT 

a(Ai = 0,r)«- — 21og(2), 
n nv 



(7) 



where we have dropped a small positive offset. The finite 
T correction to (Tmin(T = 0) is always positive, the re- 
sult ([7]) reflecting the linearly increasing density of states 
which is sampled at higher temperatures. It is impor- 
tant to note, however, that this effect does not survive in 
the presence of strong electron-electron interactions, for 
which one would obtain a length independent conductiv- 
ity at /X = 0, with a very small T dependence exhibiting 
the opposite trend. 



So far, we have taken the chemical potential in the sam- 
ple to be fixed. However, in experiments it is the charge 
density n which is controlled b y th e gate potential, rather 
than the chemical potential [l6|. As in the standard 
Fermi gas, the chemical potential is reduced upon raising 
the temperature, fi{T,n) = fiQ — {tt^ / 6)T^ / fia to lowest 
order in T <C in an infinite system. The density n and 
/io are related by Tr{hv)'^n — Hq. In general, in the ther- 
modynamic limit {W,L ^ iio/hv) the chemical potential 
/i(T, n) = Tfi satisfies 



dx ; 



2 r2 



(8) 
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withf{x,fl) = (H-exp[a; — yu])"-^. This defines the scaling 
function (1 — C (T/hv\n\^^'^Y The experimental quantity 
of interest is the conductivity as a function of T and n. 
From Eqs. (|6|8p one then easily finds the result 



h a{n, T) 



e2 L 



\ y ^ 




(9) 



My) = 2/ [c(y) + 2iog(i + . (10) 



The scaling function ^ is plotted in Fig. [TJ It has the 
asymptotics -0(0) = Vtt and ^{y I) ~ 21og(2)y, Inter- 
estingly, it attains a minimum at j/niin — 0.6932 with the 



value ip{yn 



1.5356. At large density {hvy\n\ > 



the conductivity tends to the limit 
a {\n\ > [T/hvf 



(11) 



At fixed density, the conductivity first decreases as T in- 
creases from zero (by a total of -^/tF — V'(2/min) — 13.4%) 
due to the decrease of the chemical potential. Upon 
increasing T further, a reaches a minimum at Tmin — 
J/min fi.i'lnl"'^/^ and eventually grows linearly with temper- 
ature, approaching the limiting behavior ([7]). The non- 
monotonicity in the T dependence should be observable 



for relatively low densities. 



10i°cm-2, for which T„ 



lies within the experimental temperature window, while 
being low enough for the ballistic, non-interacting ap- 
proximation to be applicable, as discussed above. Such 
a non-monotonicity has been reported in Fig. 3c of 
Ref. [l3|, whereas it was probably masked by inelastic 
scattering in Ref. At higher densities we predict a 
slight decrease of conductivity with increasing T. We em- 
phasize that the decrease of conductivity results merely 
from the decrease of fi without invoking scattering, while 
the increase of a at higher temperatures reflects the ther- 
mal sampling of the higher density of states. This is to 
be contrasted with the scenario of Ref. 15| where a non- 
monotonicity was found as well, resulting however from a 
competition between T-dependent screening and thermal 
sampling of energy dependent scattering rates. 
In order to facilitate the comparison with experimental 
data, we show in Fig.[2]the dependence of a on the carrier 
density n for various temperatures T , based on the full 
numerical evaluation of Eq. ([5]) for finite W . Apart from 
confirming the above discussed trends, we note a good 
qualitative agreement with the observations in Ref. 
at least in the regime of larger n where interactions are 
expected to be weak: the conductivity decreases with T, 
the relative decrease becoming smaller as n increases. 
Fig. [2] also illustrates the appearance of Fabry-Perot os- 
cillations (FPO) at low temperatures and finite densities, 
as predicted in Ref. 0]. This effect was not included so 
far, since we dropped the oscillatory term in ((4]) and ap- 
proximated the ribbon as infinitely wide. 
As shown in Fig. [31 the FPO's depend on the aspect ra- 
tio. Moreover, the peaks are rather easily washed out 
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FIG. 1: The conductivity as a scaling function of T /hv\n\^^'^ , 
and its asymptotics at large argument, ip[y ^ 1) = 21og(2)i/. 
At fixed n, a is non-monotonic in temperature, going through 
a minimum at Tmin = 0.693/iw|n|^''^. 
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FIG. 2: (color online) The conductivity a plotted as a func- 
tion of density n for various temperatures Q = T/Ew- The 
analytical approximation, neglecting oscillatory terms (thin 
full lines) is excellent except at low T and n. In the latter 
regime, finite size effects become visible in the form of Fabry- 
Perot resonances. [To be compared to Fig. 3c of Ref. [13] •] 



at finite temperatures. There are in fact two kinds of 
FPO's, as one can see from the inset of Fig. [3] Upon tun- 
ing the chemical potential or gate voltage, on one hand 
there are slow oscillations of periodicity A/x^ = TriS^, de- 
scribed by the subleading term in ([¥]). They originate 
from the modes which traverse the sample in relatively 
straight paths, having small transverse quantum num- 
bers n. These oscillations are washed out by thermal 
smearing when T > A/i^. However, their amplitude is 
never bigger than yj2/\\/-K « 0.14. Indeed, they can 
hardly be discerned in Fig. [5J A second type of oscil- 
lations is due to modes which enter the sample nearly 
at grazing incidence (with largest quantum numbers n). 
Those scatter back and forth between the edges many 
times when W ^ L (and are thus very sensitive to edge 
roughness). Their contribution can be obtained analyz- 
ing the transmission factors T„ as a function of chemical 
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potential fi. Fabry-Perot resonances are expected when 
there is a propagating mode at the Fermi level whose lon- 
gitudinal wavevector k is commensurate with the length 
of the sample, i.e., for 



k = = —j— , 771 = 1 , 2 , , 



(12) 



[m — does not lead to a resonance). With transverse 
wavevectors qn, such resonances occur at 



{13) 



One reads off from that the fast oscillations corre- 
sponding to fixed TO appear with roughly periodic spac- 
ings A/1/ — ttEw- In the main panel of Fig. [3l the sharp 
and soft peaks correspond to m~ 1,2, respectively. 
To estimate the width of those peaks at T = we start 
from a resonance, /i = /li™^ and fix g = g„. As increases 
by Six, k changes by Sk « S^/{tiv)^rmr. The peak due 
to the resonant transmission coefficient T„ [Eq. ^] is 
reduced to half its maximum when Sk satisfies 



{hvqnf cos^ {L{km + 5k)) 
^ihvqn)^{LSk)^/2. (14) 



This translates into a change in chemical potential of 



(m) 



%/2(m7r) 



2^L 
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The sharpest peaks correspond to to = 1. They are 
well separated form each other if 5pS^^ <C A/i/, or 
h/El > 2^/**(7rW/L)^/2, and reach an amplitude close 
to 1 at large ^. These most visible peaks start to broaden 
and decrease in amplitude when T > 5^''^\ cf., Fig. O 
Since is a rather low energy scale, and because of 

the sensitivity of these fast FPO's to edge roughness and 
inelastic scattering, their experimental observation might 
prove very challenging. 

In summary, the temperature-dependence of the conduc- 
tivity of ballistic graphene exhibits a rather unexpected 
behavior. The minimum conductivity always increases 
with temperature saturating quickly to a linear depen- 
dence on T which is a hallmark of weakly interacting 
systems. At finite density, we predict the conductivity 
to slightly decrease at low temperatures, but to increase 
again for T > T,„in ~ |n|^/2. Finally, we demonstrate two 
different types of Fabry-Perot oscillations, showing that 
these fingerprints of ballistic transport are very fragile 
with respect to finite temperatures. 
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FIG. 3: (color online) Conductivity as a function of chemical 
potential fi for various temperatures 6 = The amplitude 
of the fast Fabry-Perot resonances decays at rather small tem- 
peratures e > = V2^^^^§^ (^)^ The inset shows the 
presence of fast and slow oscillations in the deviation from 
the linear background, a = ^^^'^^^ — ^ , for a wide sample. 
The fast oscillations become more pronounced with growing 
/I, while the slow ones decrease in amplitude. 
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This is correct if the geometric capacitance per area, 
Cgeom oc 1/d {d being the distance to the gate) is 
much smaller than the electronic capacitance per area. 
Col = e'^dfi/dn = {2/-K^^'^){e^/hv)n^^^ (at T = 0). 



